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1  Foreword 

This  is  the  final  report  on  the  ARO  research  contract  DAAL03-89-G-0082- 
DAAL03-89-G-0082  issued  to  the  University  of  Maryland  Baltimore  County 
on  9/1/89  for  an  initial  period  of  3  years.  Upon  the  request  of  one  of  the  Pis 
(Rouben  Rostamian)  and  approval  by  the  program  director,  the  length  of  the 
contract  was  extended  by  a  period  of  two  years.  This  was  to  accommodate 
Rostamian ’s  leave  of  absence  from  the  university  for  a  period  of  26  months, 
where  he  served  as  the  program  director  of  the  Applied  Mathematics  Program 
at  the  National  Science  Foundation.  The  research  was  resumed  and  brought 
to  a  completion  upon  his  return  from  NSF. 

2  Project  description 

The  general  objective  of  the  research  was  to  investigate  the  propagation  of 
waves  in  stratified  elastic  media.  Of  particular  interest  were  problem  on  the 
design  of  elastic  coatings  which  due  to  their  layered  nature  deflect  the  incident 
energy  of  the  waves  and  therefore  can  affect  the  reflectivity  properties  of  solid 
surfaces  in  interesting  ways. 

3  Summary  of  the  most  important  results 

Based  on  some  preliminary  investigations,  we  made  the  following  statement 
in  the  original  proposal  submitted  to  ARO  (page  1): 

Although  no  one  coating  will  eliminate  all  reflected  or  transmitted 
waves,  we  can  optimize  a  coating  by  minimizing  the  maximum 
energy  of  a  reflection  or  transmission  coefficient  over  a  range  of 
frequencies. 

Surprisingly,  it  turned  out  that  we  can  do  much  better  than  what  we  had 
hoped  for.  The  major  result  of  our  investigation  is  that  it  is  possible  to 
eliminate  reflected  waves  of  a  prescribed  bandwidth  from  a  solid  object  com¬ 
pletely.  (See  section  4  of  this  report.) 

We  have  also  investigated  the  effects  of  air  bubbles  embedded  in  an  ab- 
sorbant  coating  layer  affixed  to  a  solid  surface.  We  have  developed  techniques 
and  implemented  algorithms  for  predicting  the  strength  per  unit  weight  of  a 
bubbly  elastic  material. 
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Figure  1:  The  frequency  response  of  a  one-layer  system. 

4  Wave  propagation  in  layered  media 

Consider  a  two  elastic  half-spaces  which  sandwich  between  them  an  elastic 
slab.  The  elastic  properties  of  the  two  half-spaces  are  supposed  to  be  known 
while  those  of  the  slab  are  unspecified.  All  three  media  are  isotropic. 

Consider  a  steady-state  harmonically  oscillating  wave  impinging  form  in¬ 
finity  perpendicular  to  the  surface  of  the  slab.  The  ratio  of  the  reflected  energy 
over  the  incident  energy  is  called  the  reflectivity  of  the  slab. 

It  is  not  difficult  to  show  that  for  any  given  frequency  w,  there  exists  a 
set  of  elastic  moduli  for  the  slab  such  that  the  reflectivity  is  zero.  Figure  1 
shows  the  frequency  response  of  such  a  slab  that  is  designed  to  annihilate 
waves  of  frequency  u  ~  20.  We  see  in  that  graph  that  when  u  deviates  from 
the  design  frequency,  the  reflectivity  rapidly  becomes  non-zero. 

To  extend  the  range  of  frequencies  over  which  the  reflectivity  remains 
small,  we  replace  the  original  slab  with  two  distinct  slabs,  each  of  which  be¬ 
ing  half  as  thick  as  the  original.  The  additional  freedom  gained  in  choosing 
two  sets  of  elastic  properties  instead  of  one,  allows  us  to  successfully  annihil¬ 
ate  two  arbitrarily  prescribed  frequencies,  as  shown  in  Figure  2.  If  we  split 
the  original  slab  into  a  larger  number  of  distinct  slabs  of  increasingly  smaller 
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< - width  of  slab - > 


Figure  4:  The  frequency  response  of  the  variable-properties  slab  shown  in 
figure  5. 


Reflectivity  versus  frequency.  Smooth  layers. 


Figure  5:  The  impedance  of  a  variable-properties  slab  that  annihilates  incid¬ 
ent  waves  of  frequencies  within  the  20-160  range. 


thickness,  we  can  impose  additional  frequencies  which  are  annihilated.  Fig¬ 
ure  3  shows  the  frequency  response  of  a  system  with  five  slabs  and  five  totally 
annihilated  frequencies. 

We  can  carry  this  process  to  its  limit  by  splitting  the  original  slab  into 
infinite  number  of  layers  of  infinitesimal  thickness.  The  elastic  properties  of 
such  a  slab  become  continuously  varying  functions  within  the  slab.  Figure  4 
shows  the  impedance  of  such  a  slab,  designed  to  annihilated  incident  waves 
in  the  20-160  frequency  range.  The  frequency  response  of  the  corresponding 
system  is  shown  in  Figure  5.  The  apparent  deviation  of  the  frequency  response 
from  zero  within  this  range  is  due  to  the  accumulation  of  numerical  errors. 

5  Optimization:  Theory  and  algorithms 

The  basis  of  numerical  computations  that  lead  to  the  optimal  impedance 
profile  of  Figure  4  lies  in  an  optimization  and  control  problem  that  can  be 
described  as  follows. 

Let  r(u)  denote  the  reflectivity  of  the  slab  corresponding  to  an  incident 
wave  of  frequency  u>.  If  we  regard  the  mechanical  properties  of  the  material  in 
slab  as  unknown,  then  the  reflectivity  depends  on  these  mechanical  properties 
as  well  as  on  u.  Since  ultimately  the  mechanical  properties  is  a  continuous 
function  of  the  depth  within  the  slab  (of  total  thickness  T),  we  regard  the 
reflectivity  as  function  r  :  [uq ,  a>2]  x  <7°[0,  T]  C  and  write  r  =  r(u>,  k).  The 
interval  [u)\,u)2}  is  an  arbitrarily  prescribed  frequency  band. 

We  wish  to  minimize  the  worst-case  reflectivity  over  the  frequency  band. 
Hence  we  arrive  at  the  min-max  problem: 

min  max  \r(uj,  «)|, 

KeKcc°[QF\u>e[u>uLJ2]1 

where  K  is  a  suitably  chosen  subset  of  (70[0,  T ]  that  incorporates  in  it  certain 
physical  constraints,  such  as  the  requirement  that  density  cannot  be  negat¬ 
ive  or  that  the  elastic  moduli  should  be  such  that  the  corresponding  elastic 
material  is  strongly  elliptic. 

The  computational  task  for  determining  the  optimal  solution  of  this  con¬ 
strained,  infinite-dimensional  min-max  problem  is  non-trivial.  The  papers 
[11],  [12],  and  [13]  concern  a  new  family  of  dual  optimization  algorithms  mo¬ 
tivated  by  the  goal  of  solving  the  min-max  problem.  By  dualizing  certain 
difficult  constraints  in  an  optimization,  much  faster  convergence  is  obtain 
than  in  competing  algorithms.  The  papers  [9]  and  [10]  analyze  the  stability 
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of  optimal  control  problems  relative  to  a  parameter.  These  stability  results 
are  applied  to  the  analysis  of  various  numerical  algorithms,  penalty  and  mul¬ 
tiplier  methods,  sequential  quadratic  programming,  and  Euler  discretizations, 
for  solving  the  optimal  design  min- max  problem. 

6  Perforated  composites 

The  presence  of  the  air  bubbles  within  an  elastic  material  reduces  the  effective 
elastic  moduli  and  enhances  movement.  That,  combined  with  viscoelastic 
effects  can  help  reduce  reflections  in  an  active  way,  by  converting  mechanical 
energy  to  heat,  as  opposed  to  simply  re-routing  the  energy  in  a  different 
direction. 

The  investigation  of  the  “bubbly”  composites  is  a  subject  of  interest  of 
its  own,  independently  of  the  issue  of  the  design  of  absorbant  coatings.  We 
have  made  some  preliminary  investigation  on  the  dependence  of  the  strength 
of  bubbly  elastic  materials  per  unit  weight  basis.  There  are  two  opposing 
effect  in  determining  the  strength  of  a  bubbly  material  per  unit  weight:  The 
introduction  of  air  bubbles  in  the  material  decreases  the  effective  density 
decreases  and  therefore  the  material  becomes  more  lightweight.  Concurrently, 
the  bubbles  reduce  the  elastic  strength  of  the  material.  It  is  not  a  priori 
obvious  how  the  strength  per  unit  mass  will  behave.  A  typical  result  is  shown 
in  figure  6.  The  three  curves  labeled  Qllll,  Q1122,  and  Q1212  depict  the 
three  effective  elastic  moduli  for  an  initially  elastic  material  which  has  been 
perforated  by  a  periodic  array  of  uniformly  distributed  circular  holes.  The 
horizontal  axis  indicates  the  relative  size  of  the  holes  compared  to  the  distance 
between  the  centers  of  adjacent  holes.  The  left  extreme  corresponds  to  holes 
that  are  small  and  far  apart  and  the  right  extreme  corresponds  to  holes  that 
are  large  and  closely  spaced  (but  do  not  overlap.)  The  gap  between  the 
curves  labeled  Qllll  and  2Q1212-(-Q1122  represents  the  degree  of  anisotropy 
introduced  do  to  oriented  structure  of  the  array  of  the  perforations. 

As  expected,  the  curves  are  all  monotone  decreasing,  indicating  the  reduc¬ 
tion  in  the  strength  of  the  solid  as  we  remove  more  and  more  material  from 
it. 

In  figure  7  we  have  normalized  the  effective  moduli  by  computing  them  per 
unit  weight  of  the  perforated  material.  Here  we  see  that  the  curves  also  mono- 
tonically  decline,  but  at  a  slower  rate.  In  particular,  Q1212  which  measures 
one  of  the  shear  moduli,  is  effectively  a  constant.  We  have  not  yet  completed 
our  investigation  of  the  perforated  media.  In  particular,  we  would  hope  to 
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Figure  6:  Effective  elastic  moduli  in  a  perforated  material  as  a  function  of 
size  of  perforations. 


Figure  7:  Effective  elastic  moduli  per  unit  mass  in  a  perforated  material  as 
a  function  of  size  of  perforations. 


Figure  8:  The  grid  for  finite-elements  computation  of  the  unit-cell  problem 
for  a  material  perforated  with  a  rectangular  array  of  round  holes. 


find  a  case  where  the  normalized  effective  moduli  reach  a  maximum  for  a 
certain  hole  size.  This  will  correspond  to  a  light-density  and  high-strength 
material. 

The  effective  moduli  for  the  perforated  materials  were  computed  by  ap¬ 
plying  the  method  of  homogenization .  In  the  core  of  the  computation  is  the 
analysis  of  the  unit  cell  problem.  Figures  8  and  9  depict  a  representative 
finite-elements  mesh  that  we  used  to  numerically  solve  the  unit  cell  problem 
and  which  in  turn  lead  to  figures6  and  7. 

7  Combustion 

The  mathematical  techniques  developed  in  this  project  for  the  analysis  of 
layered  elastic  media  have  found  an  unanticipated  application  in  the  investig¬ 
ation  of  combustion  of  inhomogeneous  solids  found  in  some  modern  artillery. 
In  a  joint  investigation  with  W.  Oberle  of  the  Weapons  Technology  Direct¬ 
orate  of  the  Aberdeen  Proving  Ground,  we  are  currently  investigating  the 
effective  chemical  characteristics  of  finely  layered  explosives. 
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Figure  9:  The  grid  for  finite-elements  computation  of  the  unit-cell  problem 
for  a  material  perforated  with  a  hexagonal  array  of  round  holes. 
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